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In systems with small spin-orbit coupling, current-induced torques on the magnetization require 
inhomogeneous magnetization textures. For large spin-orbit coupling, such torques exist even with- 
out gradients in the magnetization direction. Here, we consider current-induced torques in ferro- 
magnetic metals with both Rashba spin-orbit coupling and inhomogeneous magnetization. We first 
phenomenologically construct all torques that are allowed by the symmetries of the system, to first 
order in magnetization-direction gradients and electric field. Second, we use a Boltzmann approach 
to calculate the spin torques that arise to second order in the spin-orbit coupling. We apply our 
results to current-driven domain walls and find that the domain-wall mobility is strongly affected 
by torques that result from the interplay between spin-orbit coupling and inhomogeneity of the 
magnetization texture. 



I. INTRODUCTION 

Current-induced torques on the magnetization in con- 
ducting ferromagnets are one of the main topics of re- 
search in spintronics. In addition to being fundamen- 
tally interesting, these torques are also key to develop- 
ments in memory technology.^ Current-induced torques 
can be used to move domain walls through a ferromag- 
netic wire. When a domain wall is present the direction 
of the magnetization depends on the position in the wire. 
This spatial dependence of the magnetization gives rise 
to a mismatch between the electron spin polarization and 
local magnetization resulting in the adiabatic reactive^' 
and dissipative (also known as non-adiabatic) spin trans- 
fer torques (STTs)P^ The occurence of these two spin 
torques is well estabished but their relative magnitude, 
parametrized by the dimensionless parameter /3 which 
describes the relative strength of the dissipative torque 
with respect to the reactive one, is hard to measure^ and 
calculate.^ 

That there exist other current-induced torques related 
to spin-orbit (SO) coupling of the cariers has been pro- 
posed recentlyP^HSl In these works systems with SO 
coupling and homogeneous magnetization are consid- 
ered. Recent experiments can be interpreted using these 
current-induced torques originating from the SO cou- 
pling of the carrieraSSEHl that, unlike the adiabtic STT 
mentioned above, do not require magnetization gradi- 
ents. (Note, however, that these observations can also 
be described via the Spin-Hall effect in Pt as argued in 
Ref. [17].) For Rashba SO couping two current-induced 
torques have been found. In the experimental works a do- 
main wall is present. This implies that the description in 
terms of a homogeneous magnetization is incomplete and 
a more systematic description including both SO cou- 
pling and an inhomogeneous magnetization is called for. 

It is the purpose of this work to give such an inclu- 
sive description that incorporates both SO coupling and 
inhomogenous magnetization textures. For definiteness, 
we focus on the Rashba SO coupling. In Sec. [TT] we 
consider all current-induced torques which are allowed 



by the symmetries of the system. As the number of al- 
lowed torques is considerable, and because the symme- 
try considerations do not yield their relative magnitudes, 
we investigate these within a semi-classical Boltzmann 
description. In Sec. VI the results for the torques are 



used to calculate their effect on domain-wall dynamics. 
We find that the current-induced domain-wall velocity 
depends strongly on wall geometry. Furthermore, the 
domain-wall mobility depends strongly on the inclusion 
of torques that result from the interplay of SO coupling 
and gradients in the magnetization. 



II. SYMMETRY CONSIDERATIONS 

In this section we use symmetry considerations to ob- 
tain all allowed current-induced torques. To illustrate 
our method we begin with the adiabatic spin torques in 
the absence of SO coupling. Subsequently we investigate 
the situation with SO coupling. We use the s — d model 
since this is a convenient model to get the qualitative de- 
scription of current-induced torques. In this model the 
magnetization resides on the d-orbitals and transport is 
due to the mobile s-electrons. We investigate the system 
well below the Curie temperature, which means the mag- 
netization is represented using a unit-vector field since 
fluctuations in its magnitude are negligible. 



A. Absence of Spin-Orbit Coupling 



Within the s — d model the Hamiltonian is given by 



Usd = U {x,p) - ■ s, 



(1) 



where Hq is the Hamiltonian that descibes the motion of 
the itinerant electrons and depends on electron momen- 
tum p and postion x. We have an exchange coupling 
between the magnetic texture m(x,t) and the electron 
spin s(t) specified by the exchange splitting A. The to- 
tal Hamiltonian Ji s d is invariant under two independent 
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rotations of the spin and physical space, parameterized 
by the rotation matrices TZV and respectively. (We 
neglect the coupling between the magnetization and the 
orbit of the electrons that occurs via the Lorentz force. 
We neglect this effect for the moment because the mag- 
netic field induced by the magnetization is very small.) 
Moreover, in this description we neglect the ionic lattice. 
We explicitly have for the rotations 



K ij x j , 



P 



-s 
TZ ij p>. 



(2) 
(3) 



Note that we use the summation convention of summing 
over repeated indices. The invariance of the Hamiltonian 
implies W s d(p,x,s,m) = T-L s d{p,x 7 s,m). This means 
that these symmetries should be respected at the level 
of the equations of motion. We are interested in the 
(linear-response-) current-induced torques, hence our ex- 
pressions for the torques should be linearly dependent on 
the applied electric field E. The possible torques that are 
first order in the electric field E, which transforms under 
the action of R, should involve an inner-product with an- 
other vector that transforms under the same rotation and 
in this way creates an invariant scalar. The only other 
vector that transforms in this way for this system is the 
gradient V that acts on the magnetization. These con- 
straints lead to the two possible current-induced torques 



dm 

~dt 



oc (E ■ V)m + (3m x (E ■ V)m. 



(4) 



ST 



For a treatment of spin transfer torques that incorporates 
the symmetries of the lattice see Ref. [H] . These terms 
are frequently written in terms of the current but we 
choose to put in the electric field here as the external 
perturbation, to be consistent with the rest of this paper. 
Note the parameter (3 which is defined as the ratio of the 
dissipative and reactive spin transfer torques. 

The two torques in Eq. Q are mutually perpendicular. 
Moreover they transform differently under time reversal, 
since they differ by a factor m which is odd under time- 
reversal. This difference in behaviour under time-reversal 
symmetry implies the torques form a pair where one is 
reactive and the other is dissipative. 



B. Spin-Orbit Coupling 

?? In the presence of SO coupling the Hamiltonian 
for the spin of the s-electron couples the spin and the 
momentum of the electron. We represent SO coupling 
for spin- \ carriers via the Hamiltonian 



(5) 



where f2 contains both the exchange interaction of Eq. 
([TJ) and SO coupling, and can be seen as a position and 

momentum dependent effective exchange splitting. 

For definiteness, and motivated by experiments ] 16 * 17 ^ 
we study the simplest form of SO coupling described by 



the Rashba Hamiltonian 22 "Hr = — Xr(p x e z ) ■ s. The 
Rashba coupling together with the exchange interaction 
results in 



fl(x,p) 



-m(x) + Ar,p x e 2 



(6) 



Rashba SO coupling occurs in two-dimensional electon 
systems with inversion asymmetry along the direction 
perpendicular to the two-dimensional electron gas (which 
we choose as our z-axis). The SO coupling breaks the in- 
variance of the Hamiltonian under separate rotations of 
the spin and orbital parts of the motion. Total angu- 
lar momentum is still conserved due to the invariance of 
the Hamiltonian under combined rotations of spin and 
physical space, parameterized by TZg — IZ^ . 

The linear-response matrix L clt (m, e 2 , Vm) that de- 
scribes the current-induced torques is defined by 



m ' 



L%{m,e z ,Vm)E\ 



(7) 



where E is the electric field in the plane and e z is a 
unit vector in the z-direction. The linear-response ma- 
trix depends on this direction since inversion symmetry is 
broken along this direction. The Hamiltonian is invariant 
under parity transformations which implies that the lin- 
ear response matrix should obey —L C i t (m, — e z , — Vm) = 
L c it(m, e z , Vm). This shows that there can be torques 
on the magnetization without a gradient in the magne- 
tization. These torques tst* = Lij(m, e z )E 3 have been 
found before^ and are given by 



Tg^ oc m x (E x e z ); 



ST m x (m x (E x e z )). 



(8) 
(9) 



The spin torques are perpendicular to m because it is a 
unit vector field. Since the magnetization is embedded in 
three-dimensional space there is a two-dimensional plane 

(i) 

perpendicular to it. This means that any spin torque Tg^ 
allowed by the symmetry of the system immediately de- 
fines another torque via = mxT^. These pairs dif- 
fer a factor m which changes its sign under time-reversal, 
hence the two torques form a reactive-dissipative pair, 
like the STTs in Eq. @. In the following we will show 
only one of the pair. All terms to first order in the gradi- 
ent of the magnetization that do not involve e z , are given 
by 



_(2) 
X ST 


oc 


(E ■ V)m; 


(10) 


T (3) 
X ST 


oc 


((m x E) ■ V)m; 


(11) 


T (4) 
X ST 


oc 


(m ■ E)(m ■ V)m; 


(12) 


T (5) 
X ST 


oc 


E a (m x V)m a ; 


(13) 


_(6) 

x st 


oc 


(m x E) a (m x V)m a ; 


(14) 


_(7) 

x st 


oc 


rn x E(V ■ m); 


(15) 


T (8) 

x st 


oc 


(m x E)m ■ (V x m), 


(16) 



In the first line the familiar STT^S describing the 
current-induced torque in systems with inhomogeneous 
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magnetization is obtained. Together with the dissipative 
g'P'jEHSI that is associated with it (rf^) those torques de- 
scribe the weak SO coupling situation. In the second line 
we find a STT due to a Hall current. The other torques 
do not have a straightforward physical interpretation. 

Up to this point we have explicitly given the torques 
to first order in either e z or V. There are more torques 
that involve an even number of e z 's and are first order 
in V. We will not list them because the list will be to 
long to be illuminating. We proceed by actually calcu- 
lating the torques in the next section. The reason we 
do this is twofold. First, having demonstrated the exis- 
tence of many spin torques due to the combined effects of 
SO coupling and magnetization gradients, we now explic- 
itly calculate which torques occur within a semi-classical 
approach to the Rashba model. The second reason is 
to give an estimate of the relative magnitude of the vari- 
ous current-induced torques which cannot be found using 
symmetry arguments. 



III. SEMI-CLASSICAL FRAMEWORK 

In order to investigate microscopically which current- 
induced torques appear for the textured Rashba ferro- 
magnet we use a semi-classical approach. This approach 
has proved its merit in the description of the anomalous 
Hall effect PS We describe the system by the Hamiltonian 



H = 



2m„ 



fl(x,p) ■ s + E MM [m] 



(17) 



where fl(x,p) is the effective Zeeman field, given in Eq. 
([6]), that incorporates the Rashba SO coupling and the 
exchange coupling, and Emm [m] is the micromagnetic 
energy functional for the magnetization. Furthermore, 
m e is the effective mass of the electron. The equation of 
motion for the spin degree of freedom is written as 

ds 1 _ a . . 
— - = —s x 12 — — s x (s x Si), 
at h h 

where we introduced a damping term proportional to a 
that describes relaxation of the spin into the direction of 
the effective Zeeman field. The spin dynamics is much 
faster than the motion of the electrons such that we can 
solve the above equation of motion up to first order in 
time derivatives of fl. We obtain the following solutions 

a h dfl - ha dfl , „. 
s s = sfl + s . — xfl -. — , (18 

Vfi^fi dt dt 

where s = ±1 describe the majority(s = 
+) /minority (s = — ) electrons, and fl = 0/|0|. 
The first term describes the adiabatic following of the 
effective magnetization texture by the electron-spins. 
The other terms describe the slight mismatch of the 
spins with the effective magnetization. We find the 
dynamics of the itinerant electrons by inserting the first 



order solutions of the spin degree of freedom, given in 
Eq. (18), into the Hamilton equations of motion for the 



electrons. We obtain 



de^ 
dp 1 



sh 



dfl dfl 



dp 1 



dt 



de s 
dx l 



dfl dfl 

d~x l * ~dt 



fi + ah 



fl 



dfl dfl 

dp 1 dt 

_dfl dfl 

dx l dt 



\e\E, 



where e s = p 2 /2m e — s\fl\ is the dispersion for the 
majority (s = +) /minority (s = — ) electrons. Note that 
we added an electric field to induce a transport current. 
The total time derivatives on ft should be understood as 



dfl 

~dt 



dfl 

: d~x l 



' Ps 



dfl 

dp 1 ' 



Now that we have this semi-classical description of the 
single particle dynamics we calculate the spin-torques us- 
ing the Boltzmann equation for the distribution function 
f s (x,p,t) of the particles, which, in the relaxation-time 
(r r ) approximation, is given by, 



dt 



fs(x,p,t) 



f.(x,p,t)-f™(e B ) 



(19) 



where / FD (e) = (1 + e^ e ) _1 is the Fermi-Dirac distri- 
bution function. The relaxation-time approximation is 
the simplest description of the Boltzmann collision in- 
tegral. We make the relaxation-time approximation here 
for convenience. A detailed study of the collision-integral 
in the presence of strong SO coupling is beyond the scope 
of this work. We refer to the work by Pesin and Mac- 
Donald in Rcf. pjj] for more details on the situation of 
homogeneous magnetization. The left-hand side in Eq. 
( 19 1 should be read as 



df s = df s (x,p,t) ^ | df s (x,p,t) . 
dt dm dp 



df s (x,p,t) 
dx 



The equation of motion for the magnetization is the 
Landau-Lifshitz-Gilbert (LLG) equation 



dm dm 
— - = -7m x ifcff + a G m x — — 
dt dt 



T~sd, 



(20) 



where 7 is the gyromagnetic ratio and the torques due 
to the s — d coupling r s d = A/(2fi)m x (s) contain the 
spin torques of interest and a renormalization of the pa- 
rameters in the LLG^ equation which we discuss in this 
section. The current-induced torques are proportional to 
the electric field and will be given in Sec. |V| The renor- 
malized LLG equation we obtain is given by 



(1 - n) 



dm 

~dt 



-jm x H' cS + a' G m x 



dm 

~dT 



TST, (21) 



where Tst contains all terms of t s( j proportional to the 
electric field and H' oS is defined as the effective magnetic 
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field acting on the magnetization which acquires an ad- 
ditional term from the coupling to the electrons 

^ ff = ^T +a ^ i(1+ A^ )(m ' ez)e - (22) 



and the rcnormalized quantities in Eq. (21 1 are given by 



m e a ( A m e \ ep 2 

— i^r — 4 t m y 

irh 2 \ 2 A 



where a is the lattice constant. The additional term in 



Eq. (22 1 is an anisotropic damping term which for the 



typical parameters (see Table |H|) we use in the calculation 
of the domain wall dynamics is negligible, moreover these 
parameters also imply 77 <C 1. Furthermore we obtain 
that the observed Gilbert damping constant is given by 



a 2 m e 
Irf 2 



^-^(1-^(1 + 4^)) 



(23) 

Note that ag phenomenologically describes the damping 
of the magnetzation due to interactions other than the 
s — d coupling, such as relaxation due to magnon-phonon 
interactions. Before we calculate the spin torques within 
this semi-classical framework we determine the current as 
a function of electric field within this simple model. We 
need this later on to express the spin torques in terms of 
the current. 



where Go = 2\e\ 2 /h is the quantum of conductance. The 
second contribution 

^ = | 2 ^ 2 m z +Ar r (m.(Vxm)-«|(V.m))j f «^ 

,( 26 ) 

which is the anomalous Hall response generalized to inho- 
mogeneous magnetization. The last contribution to the 
conductivity is 



g AMR _ x2 ( 4e F T, 



( 4eF T r \ iaz jbz^ab / 7\ 
A 2 " ~ T ) ' ( 27 > 



which depends on the relative orientation of the electric 
field and the magnetization and hence corresponds to 
anisotropic magnetoresistance. We also define the cur- 
rent polarization via 



for later reference. 



E 



d 2 P 

iJ (27^ 



f(e s )sx s , 



V. CURRENT-INDUCED TORQUES 

In this section we give the current-induced torques for 
the Rashba model, introduced in Sec. |III| The current- 
induced torques can be calculated from the current- 
induced spin density. They are given by 



IV. CONDUCTIVITY 

In this section we give the conductivity for the Rashba 
system. Note that the conductivity we find here is only 
correct within this simple s — d description. We need 
the conductivity in order to interpret the current-induced 
torques in the next section. The conductivity cry is de- 
fined as j l c = a % i E l , where j c is the charge-current den- 
sity. We calculate the conductivity up to first order in 
the gradient of the magnetization and up to second or- 
der in the SO coupling strength. The expression for the 
charge-current density is given by 



E 

s=± 



(2tt/i) ; 



f s (x,p,t)x s 



(24) 



Using the relaxation-time approximation described in the 
previous section we find that the conductivity has three 
contributions a = <to + ^ah + camr corresponding to the 
diagonal, anomalous Hall effect and anisotropic magne- 
toresistance, respectively. The diagonal conductivity is 
given by 



Go 



e F T r 



4m e X 2 eF 



8 l \ 



(25) 



1"ST 



2h 



(s) 



2h 



-m x 



E 

s=± 



d 2 p 



f s (x,p,t)s s , (28) 



where tst is the sum of all the separate spin torques we 
list below. We evaluate the integral in Eq. |28| ) up to first 
order in the damping parameter a and gradient of the 
magnetization and up to second order in the spin-orbit 
coupling strength A. Note that we only include terms 
linear in the electric field, and that we take dm/dt = 0. 
Taking into account this time-dependence gives rise to 
renormalization of damping and gyromagnetic ratio that 
we already discussed in the previous section. 

In agreement with our phenomenlogical arguments [see 
Eqs. (|8jrefeq:tstlp)], we obtain two spin torques that are 
zeroth order in the gradient of the magnetization which 
are given by 



-(i) 



e\m e \a 2 ( Ar r 2ep 



nh 2 \ 2h 
,2 



a— ] (E x e z ) x m; (29) 



n 1 \ \e\m e \a ,,_ . . 
= t 2 (( x e - x m 



(30) 



These homo geneo us SO induced spin torques where de- 
rived before.Q^EH j n case a = 0, we agree with Man- 



chon and Zhang 12 and with Kim et. al™ about the ratio 
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between the reactive and dissipative torques. In addi- 
tion, for a / 0, we find another contribution to these 
torques coming from spin relaxation. Note that the two 
torques given above form a perpendicular pair, one dis- 
sipative one reactive. In what follows we will group the 
torques into these pairs when both reactive and dissipa- 
tive torques emerge to second order in SO coupling. 

All other torques are first order in the gradient of the 
magnetization. The first two torques we find are given 
by 



_STT 



„STT_L 



■Kh 2 V 2 r r A A 
(E ■ V)m; (31) 

^ 2 (j+|l(3 + 4^))) 
(32) 



\e\ar r a / m e . 
m x (E ■ V)m, 



which are the well known STTs that also occur in sys- 
tems with negligible SO coupling, see Eq. Q, and are 
due to the spin-polarized current in the direction of the 
electric field. The ratio of these two torques defines the 
j3 parameter. We find that 



2a 

" a" 



A 



r (30 + 8m 2 ) 



(33) 



In the previous section we showed that the current can 
be decomposed into three components. Several of the 
torques we find can be interpreted as the ordinary spin 
transfer torques [Eq. Q] with current response modified 
due to the SO coupling. First, we have the torques 



4\e\m e e F \ 2 a 2 



irA 2 h 
2\e\m e a\ 2 a 2 



(m • e z )((E x e z ) ■ V)m; 



m ■ e z )m x ((E x e z ) ■ V)m, 



which are due to the anomalous Hall current, 



a l p^E 3 , and are given to first order in gradients as 



.AH 



.AHl 



= V{jab. ■ V)m; 

= P~^— m x Wah • V)m. 



(34) 
(35) 



Two torques can be interpreted to be a generalization of 
the STTs coming from the anisotropic magnetoresistance 
response given by Eq.(27). These torques are 



.AMR 



2|e|m P A 2 a 2 



nAK 

.amri _ \e\m e X 2 aT r a 2 

■nh 2 



(a - 24^) ({E x e z ) • m)((m x e z ) ■ V)m; 



(5 + 16-£r)((E x e z ) ■ m)m x ((m x e z ) ■ V)m. 
A z 



(36) 
(37) 



The next torque, given by 
\e\m„X 2 aT r a 2 



.Hall 



2ttH 2 



(1- 



-A^)((Exm)-V)m, (38) 



has the symmetry of a STT due to a normal Hall 
resonse. This is not the normal Hall response because 
it is quadratic in the SO coupling parameter. In our de- 



scription we did not include the normal Hall response of 
the system, due to the smallness of the effect. 

The torques obtained up to this point could be inter- 
preted as the known SO coupling induced spin torques 
for Eqs. ([29] |30j) and the STTs [in Eqs. Q-g] with 
a current response that is modified due to SO coupling. 
Now we will list the torques that cannot be interpreted 
as known current-induced torques. We have the pairs 



a 2\e\m e \ 2 a 2 fepT, 
■nAh 

aJ _ 4\e\m e e F \ 2 a 2 

= 7rA 2 /l 

.a|e|m e A 2 a 2 r r 
wh h 
2\e\m e X 2 a 2 



= -2 



-b_L 



— a ) ^ m x ^ x ez ^ am x ( ez x ^) mQ ' 

(m x (E x e z )) a m x (m x (e z x V))m a , and 
(E x e z ) a m x (e z x V)m a ; 



irAh 



( a ~~ ~Ji) ^ E X e2 ^ m x ( m x ( e * x V))m Q . 



(39) 
(40) 
(41) 
(42) 



We also have four torques that do not form reactive- dissipative pairs we list them below 

\e\m e €FX 2 T r a 2 



t c = 4^ 



nAh 2 
\e\m e €FX 2 T r a 2 



tAH 2 



(m ■ Vm • e z )m x (E x e z ); (43) 
(E x e z ) a (m ■ V)m Q (m x e z ); 



(44) 



G 



\e\m e \ 2 a.T r a 2 

3\e\m e X 2 aT r a 2 
2nh 2 



1 



A-^)E a {m x V)m a ; 



(45) 



(1 



4^K(mxe 2 )(£;.VK.(46) 



Note that the above torques are of second order in e z , and 
have therefore not been explicitly written down in Sec. 
??. The current-induced spin torques in this section are 
the central result of this paper. From the list of torques 
we presented here it is clear that the interplay of SO 
coupling and an inhomogeneous magnetization gives rise 
to many spin torques. In the next section we consider 
their effect on current-induced domain-wall motion. 



VI. DOMAIN- WALL MOTION 

In this section we investigate the effect the spin torques 
have on current-induced domain- wall dynamics. We 
study the domain-wall dynamics by employing the one- 
dimensional rigid domain-wall model. Within this model 
the dynamics is captured by the collective coordinates 
of the wall which are its position r dw and central an- 
gle <Pdw- We study three different realizations of domain 
walls summarized in Table [IJ Due to the SO coupling 
the current-driven motion of the three walls differs. In 
order to arrive at the equations of motion for the col- 
lective coordinates we describe the direction of the mag- 
netization m — (cosvJdwSin^dwjSin^dwSin^dw^os^dw) 
using two angles 9 dw and ip dw . We use 8 dw (x,t) — 
2arctan [exp(x — r dw )/Adw] and f dw (t), where A dw = 
\J J IK is the domain-wall width in terms of the exchange 
stiffness J and the easy axis anisotropy K. The direc- 
tion of the electric field is specified by th e angl e 4>e with 
the x-axis in the x-y plane. The knowrPSEU equations 
of motion for the collective coordinates rdw and <p dw are 
augmented by terms obtained from the current-induced 
torques of the previous section. In the calculations we 
make use of the parameter values as given in Table [IT] 
These parameters ar typical for metallic ferromagnets, 
and the value of the spin-orbit coupling is taken from 
Ref.[ ll6l ]. Furthermore, we give the results as a function 
of the critical field E c and velocity v c for the case without 
SO coupling, which are defined as^ 



K± 



E r = 



(47) 
(48) 



TABLE I. Magnetic anisotropy configuration and the corre- 
sponding domain wall structures. 

Easy Axis (K) Hard Axis (K±) 



Bloch(z) 

Neel(x) 

Bloch(y) 



z 

X 

y 



TABLE II. Parameters used in domain wall motion calcula- 
tions. 



CF 


= 1 eV 


A 


= 0.1 eV 


m\ 2 


= 9 meV 


a'a 


= 0.05 


a 


= 0.05 


T 


= 30 fs 




= 10 nm 


a 


= 0.3 nm 



where the spin velocity in absence of SO coupling is de- 
fined as 



2tt/j 2 



which is the zero S) coupling (A — > 0) limit of Eq. (31 1. In 



Eq. (21) we introduced the renormalized Gilbert damp- 
ing parameter a' G which is the Gilbert damping para- 
mater that will be measured in experiments. We expect 
that the Gilbert damping ac for the magnetization and 
the damping a for the itinerant spins are of the same or- 
der of magnitude. In the Appendix we give the equations 
of motion for the Neel(x) and Bloch(y) wall configura- 
tions. Here we explicitly address the Bloch(z) wall. 

The equations of motion for the collective coordi- 
nates are obtained by inserting the Bloch(z) domain- wall 
ansatz, as given above, into the equation of motion for the 
magnetization, see Eq. (21 ). To get the equations of mo- 
tion we take the inner-product with 5ra.Bioch(z) / '&*dw! for 
one equation of motion and similar for #J7iBioch(z)/^¥'dw 
Subsequently we integrate those two equations over all 
space. The two equations of motion we obtain in this 
way are given below 



r dw , K x . n f n M) ■ 6ra + 4rC 



, a G (p dw — sin 2ip dw = -t { > H — cos cp dw - — — sin ip dw E cos(0 B - (p dw ) 

Adw rl \Z OAdw oAdw / 

2r b_L _ 4r AMR r STT 

H ttt Esm((f> E ~ ip dw ) sin</?dw t £cos0_e; (49) 

OAdw ^dw 
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<^dw + «G7^ = ( - T COS</?dw ) EcOs((f>E - Ifdw) 

-"Mw \ 1 ^dw / 

4r AMR_L + 2r b 2r f + 3r STT - L 
H Esm(§ E - v? dw ) sin^ dw H — Ecos<j) E , (50) 





0. 0.02 0.04 

-E/E c 

FIG. 1. Average velocity of a Bloch(z) wall as a function of the 
applied field. The dashed (blue) line is the situation without 
spin-orbit coupling, the dotted (black) line shows the results 
with only the homogeneous SO torques, i.e. t' 1 ' and r' ix ' 
added. The solid line (red) shows the result of the solution 
of the equations of motion including all spin torques. The 
parameters used to obtain these results are given in Table [ill 




0. 0.02 0.04 

-E/E c 



FIG. 2. Average velocity of a Neel(x)-wall as a function of 
applied electrix field in the ^-direction. Lines are as in Fig[l] 
The equations of motion can be found in the Appendix. 

The scalars rW are defined as the prefactors in front of 
the vector quantities of the torques in section [V] In Fig. 
[l]we show the average Bloch(z)-wall velocity as a func- 
tion of the applied electric field in the x-direction. The 
boundary conditions for the current through the ferro- 
magnet are such that only a current in the x-direction 
is present. In the figures we took <Pe — 0, since the 
off-diagonal contributions in the conductivity give rise to 



0. 0.02 0.04 

-E/E c 

FIG. 3. Average velocity of a Bloch(y)-wall as a function of 
applied electrix field in the x-direction. Lines are as in Fig[T] 
The equations of motion can be found in the Appendix. 

a small (< 1% of the external field for the parameters 
used) voltage gradient in the y-direction. The average 
domain- wall velocity is defined as v^w — (^dw): where 
the brackets denote a long-time average. 

From the results in Figs. [T]-[3]we see that the inclusion 
of spin torques due to the combined effect of do cou- 
pling and magnetization changes the domain-wall mo- 
bility /idw = dvdw/dE completely as compared to the 
situation without these torques. 

In Figs. [2j [3] we show the results for the Neel(x) and 
Bloch(y) walls respectively. It is clear that also in this 
case the additional torques induce qualitatively different 
behaviour of the domain wall motion compared to the 
situation with only the torques induced by SO coupling 
for homogeneous magnetization. 



VII. DISCUSSION 

In this paper we considered current-induced torques 
in systems that have SO coupling and a textured mag- 
netization. The effects of these torques on domain-wall 
motion have been investigated. We have shown that the 
effects of the interplay between the SO coupling and the 
gradients in the magnetization are qualitatively impor- 
tant for domain-wall dynamics. In particular, we showed 
that the inclusions of all torques typically changes the 
domain-wall mobility as compared to including only the 
spin transfer torques that occur at weak spin-orbit cou- 
pling and/or the homogeneous spin torques due to SO 
coupling. The results of this work may be used to dis- 
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criminate between Rashba SO coupling and injection of 
a spin current via the spin Hall effect, because the latter 
will only show the homogeneous current-induced torques. 

In this paper we considered Rashba SO coupling. Our 
results can be generalized straightforwardly to linear 
Dresselhaus SO couplingpS which is linear in momentum 
too. For linear Dresselhaus coupling the dispersion of the 
carriers is the same as for Rashba coupling. The effective 
magnetization for the Dresselhaus SO coupling is given 
by ttv(x,p) — Am/2 + \£,(—p x ,p y ,0) T . This means 
P x e z -> ( 

~PxiPy>ty T when we go from the Rashba to 
the Dresselhaus coupling. The current-induced torques 
we found in Sec. |V| involve factors v x e z , where v is a 
vector. For clarity we consider r' 1 ) cx (Exe z )xm (given 
in Eq. ( 29 ) ) , for the Dresselhaus system the torque would 
be in the direction (E x e z ) x m — >• (—E x ,E y ,0) x m. 
In this way we obtain the results for the textured Dres- 
selhaus ferromagnet. The results for combined Rashba- 
Dresselhaus SO coupling are less straightforward to ob- 
tain since the dispersion of the carriers changes. 

Another obvious place to look for the appearance of 
addiational torques due to SO coupling would be in di- 
lute magnetic semiconductor systems, where the effective 
Hamiltonian for the carriers also has strong SO coupling. 
In Ref. [26] spin torques for the dilute limit are calculated 
for this system. In that work one of the current-induced 
torques is interpreted as an anisotropic dissipative STT. 
This anisotropic torque can as well be interpreted as the 
torque given by Eq. (151. It would be very interesting to 



see which other torques would appear in those systems 
because the allowed spin torques would be listable since 
there is no symmetry breaking in the z-direction. 



The reciprocal physical mechanism associated with 
current-induced torques are currents driven by non- 
equilibrium magnetization dynamics, often referred to as 
spin- motive forces. We obtain these using the Onsager 
reciprocal relations.^ We do this via the linear response 
matrix 



7 ^ f=ijk 



LH t (m,e z ,Vm) 



H 



off 



\fcj \Ls mf (m,e z Vm) <t« (m,e„Vm) J \ E' J 

where L^ t (m,e z ,Vm) is the (3 x 3) matrix that gives 
the current induced torques as defined in Eq. ^ and 
L % JLf (m, e z , Vm) gives the spin motive forces. These two 
matrices are related via Onsager reciprocity which yields 



e z ,Vm) = L J sm{ (-m,e 



-Vm) 



In future work we intend to explore these spin motive 
forces in more detail. Another interesting direction for 
future research is the inclusion of thermal gradients and 
heat currents. 
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Appendix A: Different domain-wall configurations 

In this appendix we give the equations of motion for the 
Neel(x) and Bloch(y) domain-wall configurations. The 
magnetic anisotropy configuration corresponding to these 
different walls is given in Table [I] 



1. Neel(x) Wall 

The Neel(x) wall is parameterized as m = 
(cos(6*(x, t)), cos 4>(x, t) sm.{6(x, t)), sin(j)(x, t) sin(0(x, t))) T . 
The equations of motion are obtained as explained in 



Sec. VI of the main text. The equations of motion for 



the collective coordinates are given by 



r dw , . Kx ... /ili 

a G tp dw = —r- sm2(y5 dw — I t v ; 

Adw " 



7T 



(1) 



1 



3A 



dw 



— — r c cos 2 (p dw + -J- (r AH - r b + r Ha11 ) sin <p dw ) E sin < 

•J-^dw 4Adw 



r_ 3r AMR + r bX +r d ) 



1 



3A, 



d w 



-(3r a + r c )sin^ (^ dw Ecos(j) E 



<^dw + O.Q 

^dw 




Qr e - AdwT^^M sin</?dv 



((2r AMR + r b ) cos 2 <p dw - (-3r a± + r f ) sin 2 ^ dw ) + T STT± }Ei r, t: 



-d) 



48A dv 



(4r AH1 - 4r b± - r c - r d ) sin<p dw I Esmy E . 



2. Bloch(y) Wall 



For the Bloch wall the magnetization is parameterized as m = (cos <ft(x, t) sin(9(x 7 1)), cos(9(x, t)), sm(f>(x, t) sin(9(x 7 t))) T . 
The equations of motion are 



^dw , . Kj_ . 

a G (p dw — sin 2ip dw 

Adw n 



Vdw \ 



AMR 



3Adw 
2 

3Adw 4A 



2r b - 3t stt - l - r f cos 2 ip dw + 3A dw T (1) ) E cos <j> E 



dw 



AH a_L , _Hall 

r — r 



+ T HaU + 2AdwT (1) )cos</?dw] Esin(j> E 



l ?*dw 

Vdw + OLq 

Adw 



3A d 



— f T AMR± + 2r b + 3AdwT (l) _ 3t STT± _ T f CQS 2 \ E CQS ^ 



+ 32AT + 8rAH± + r " + ^ ~ 16Adwr(1±) ~ (T C + T d ) cos2 <Aiw) cosip dw Esm0 E . 



